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We numerically investigate the spectral statistics of pseudo-energies for the unitary network 
operator U of the Chalker-Coddington network. The shape of the level spacing distribution as 
well the scaling of its moments is compared to known results for quantum Hall systems. We 
also discuss the influence of multifractality on the tail of the spacing distribution. 
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§1. Introduction 

The Chalker-Coddington network model, ^ representing systems of independent two-dimensional 
(2D) electrons in a strong perpendicular magnetic field and smooth disorder potential, is a conve- 
nient tool to investigate the quantum Hall universality class. 1_3 ) It has been used to determine var- 
ious quantities characterizing the localization-delocalization (LD) transition taking place between 
the quantized plateaus of the Hall conductance. 4 ) The first approach was to use the transfer matrix 
method to numerically determine the characteristic exponent v ~ 2.35 for the divergence of the 
localization length £ oc \E — E C \~ V when the energy of the electron E approaches the critical energy 
E c . l ^"> At the critical energy itself the network model was used to calculate critical wave functions. 
Those were subjected to a multifractal analysis and the critical exponent an ~ 2.27 describing 
the scaling of typical value of their squared amplitude with the system size, expd^d 2 } oc L _Q °, 
was determined. 2 ^ The model can also very conveniently be used to simulate the diffusion of wave 
packets using a unitary time evolution operator U. By doing so it was possible to determine the 
scaling behavior of the local density of states. 6,7 ) 

Recently, a new approach in the investigation of network models was suggested. It was argued 
that network models can be used to determine the statistics of the eigenvalue spectrum of their 
underlying system by investigating the spectrum of the unitary network operator U. 3 ^ So far 
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the numerical investigation of the pseudo-energy statistics has been concentrating on the number 
variance £2(^0 = {(n — (n)) 2 ) of an energy interval containing on average N = (n) levels and 
the compressibility x = lirn7v->oo lmi.L->-oo ^ was f° un d that the prediction by Chalker et al 8 ^ 
that x = (d — D2)/{2d), where D2 is the correlation dimension, can be confirmed. Furthermore, 
it was seen that the level spacing distribution function P(s), denoting the probability to find two 
consecutive levels separated by an energy e = sA, A being the average level spacing, deviates in 
its tail from the Wigner surmise. 3 ) This was expected at the critical energy 9 n ) and was also seen 
in other numerical investigations. 12 ~ 14 ) In this paper we are going to take a closer look at the 
level spacing distribution P(s). We will show that by using the spectrum of the unitary network 
operator U all the known results regarding -P(s) can be reproduced. Furthermore, we will introduce 
an improved numerical method to determine that spectrum by reducing the size of the matrix that 
has to be diagonalized. 

§2. The Network Operator and its quasi energy spectrum 

The Chalker-Coddington network is a 2D regular lattice whose nodes consist of elementary scat- 
terers. The scatterers are connected by unidirectional channels. Each scatterer is a 2 x 2 matrix 
S k mapping 2 incoming channel amplitudes onto 2 outgoing channels. In order to simulate ran- 
dom channel lengths, a random kinetic phase is added to each channel amplitude. Vertical links 
are always scattered to horizontal links and vice versa. The scattering probability to the left and 
the right are determined by a transmission amplitude T and a reflection amplitude R = 1 — T, 
respectively. This chiral structure of the network captures the basic features of the motion of 2D 
electrons in a strong magnetic field and smooth random potential, i.e. the potential correlation 
length d is much larger than the magnetic length I. In this picture the network links represent the 
equipotential contours of the random potential and its nodes are the saddle points. The scatter- 
ing amplitude is determined by the electron energy E and the energy Uk of the respective saddle 
point, T = (1 + exp((£ l — u^)/ 'Et))" 1 , where Et is the characteristic tunneling energy of the poten- 
tial. 2 ' 15 ) The saddle point energies can be either fixed = for all k) or randomly distributed 
(ui £ [— W, W]), in any case an LD phase transition occurs when E approaches the critical energy 
E c = 0. 1 ) In the following we assume that all saddle points have energy u& = , i.e. W = 0, 
corresponding to the original Chalker-Coddington model. This is done to avoid effects caused by 
classical percolation. We will discuss the case W > and the corresponding percolation effects in 
a later publication. 

A state or wave function of the network is defined by a normalized vector ^ = (ipi, . . . ,ip n ) 
whose elements ipi denote the complex amplitudes on the n = 21? network channels, where L is 
the system size. Its time evolution is determined by the scattering matrices at the saddle points 
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which can be incorporated into a single unitary network operator U : 

2L 2 

Mt + i) = Y,u ij i> j (t), 
j 

where the Uy are either appropiately arranged matrix elements of the scattering matrices S k or 
zero. This results in U having only two non-zero elements per matrix column and being energy 
dependent through the energy dependence of the scattering amplitudes. An energy eigenstate VP of 
the underlying electronic system is a solution of the stationarity equation: 

U(E)V = V, (1) 

which corresponds to an eigenstate of U(E) with eigenvalue 1. Such an eigenstate will occur only 
at discrete values of E = E n which form the energy spectrum of the modeled system. In contrast 
to these energy eigenvalues, the eigenvalue spectrum of U(E), defined by the equation 

U(E)V a = e^W**, (2) 

for a fixed value of E is much easier to determine. 2 ) As was pointed out in Ref. 3, the eigenphases 
uj a (E) show the same statistics as eigenenergies E n close to E. They can be considered an excitation 
spectrum at energy E and we will refer to them as pseudo-energies. 

The fact that we can choose the energy E at which we want to investigate the level statistics 
makes it possible to look directly at the critical energy E c without having to worry about leaving 
the critical region. This means that all the critical pseudo-energies, i.e. the eigenvalues of U(E C ), 
can be used for our statistics, which is an enormous advantage over other methods that can only 
use part of the eigenspectrum of their Hamiltonians. In order to further improve the efficiency of 
the numerics, we also use the fact that for any eigenphase u> a the phase oo a +ir is also an eigenphase. 
This pseudo-degeneracy is caused by the chiral structure of the Chalker-Coddington model and the 
resulting fact that each time step maps horizontal to vertical links only and vice versa. If you look 
at U 2 this property results in the formation of two orthogonal [7 2 -invariant subspaces that have 
the same eigenvalues. This U 2 degeneracy is reflected in the pseudo-degeneracy in U, which forced 
us to disregard half the eigenvalues in the past. Now we will use it to reduce the size of our matrix. 
When diagonalizing the matrix U 2 , instead of U, and calculating the degenerate eigenphases 2uj a , 
we can take out all matrix elements corresponding to vertical (or horizontal) links and thereby 
reduce the matrix to a quarter of its size. Diagonalizing this stripped matrix and dividing all its 
eigenphases by 2 gives us the upper semi-circle of all eigenvalues of U, which is the maximum yield 
possible, within a much shorter time. Now let us take a look at the results. 

§3. The level spacing distribution P(s) and scaling 

The level spacing distribution in the localized regime, i.e. at energies \E\ 3> E c = is Poissonian: 
Pp(s) = exp(— s). The closer we get to the critical energy the more it will change and become closer 
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to the GUE distribution of the metallic phase known from random matrix theory, 16 ) Pgue{ s ) = 
||s 2 exp ^— |s 2 ^ (see Fig. 1). Since we do not have a metallic phase in the quantum Hall regime, 
the GUE distribution will not be reached. At the critical energy the distribution is neither GUE 
nor Poisson, because the critical wave functions are neither smeared out homogeneously nor are 
they localized on a small area. Instead, they form a self similar measure fluctuating strongly on 
all length scales, which is best described in terms of multifractality. 17 ) It has been found that at 
criticality the P(s) distribution is best overall fitted by the following formula: 14 ) 

P(s) = AsPexp(-Bs a ), (3) 

where A and B are determined by the normalization conditions (1) = 1 and (s) = 1. This shape 
was proposed by analytical theories that predicted that a = 1 + l/{dv) ss 1.21. 10 ' 11 ) We will see 
that this prediction is not supported by our data, instead an earlier prediction 18 ) for the tail region 
of P(s) that makes a connection to the compressibility x can be verified, i.e. P(s) cx exp(— s/(2x)). 

In our investigation we determined pseudo-energy spectra for system sizes L ranging from 16 to 
50, i.e. L 2 saddle points and eigenvalues per system. The number of realizations per system size 
was set to get from 4 x 10 4 up to 5 x 10 6 eigenvalues. This resulted in enough statistics to get good 
results for P(s) for values of s up to about s = 4.0. In Fig. 1 P(s) is shown at different energies 
ranging from E = E c = to E = 3. The change from the critical distribution at E = 0, which 
does not fit the GUE behavior, to the Poissonian type, which is almost reached at E = 3, can be 
seen very clearly. This change is expected to show one-parameter scaling for any function Z(E, L) 
that describes the shape of P(s): 

Z(E,L) = f(L/0, i = Ho\E-E c \-». (4) 

As a scaling variable we consider several quantities, first 19 ) 

roo 1 rco 

J (L,E)= Q ( s )ds = - s 2 P(s)ds, (5) 
JO 2. Jo 

where Qo(s) is the probability that an energy interval of width s contains no energy eigenvalue. 
Note that Jo is half the lowest order momentum of P(s) that is nontrivial and is also connected 
to q = (s 2 )^ 1 that describes the 'peakedness' of a distribution. 20 ' 21 ) Furthermore, we invoke the 
entropic moment of P(s), — (sln(s)) and calculate, apart from — ln(q), another generalized Renyi- 
entropy 5 str = -(sln(s)) + ln(s 2 ). 20 " 22 ) 

All of these quantities are numerically easily obtainable and they take all calculated neighboring 
level spacings into account. Since all of them show nice scaling behavior we present only that of 
Jo- The inset in Fig. 2 shows Jo as a function of E for different system sizes. As expected, all the 
curves meet in a single size independent point Jq = 0.61 ± 0.01 at the critical energy E = E c = 0. 
If we assume the validity of Eq. (4) we can rescale L so that the data points of Jo as a function of 
L/£ for consecutive energies E overlap with each other. The result of such a rescaling is shown in 
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Fig. 2. In this case we fit a fourth degree polynomial in ln(L/£ ) + v\nE to a log-log plot of our 
data and thereby determined v = 2.1 ± 0.3. This is consistent with other results for i> obtained by 
various previous studies. 23 ~ 27 ) 

Now, let us proceed with the shape analysis of the P(s) function using the other two quantities 
q and S str . As it has been shown in a number of cases 20 ' 22 ) these parameters apart from showing 
scaling and fixed point in the critical regime, as a bonus enable us to determine the one-parameter 
family of P(s) functions describing the transition from e.g. GUE-like to Poissonian behavior. To 
this end we apply the same linear rescaling as was done in Ref 21 and introduce Q and S that span 
the whole [0, 1] interval with (1) corresponding to the Poissonian (GUE-like) limits, respectively. 
In Fig. 3 the parametric dependence of S vs Q is shown. We clearly see, that all the data fall on a 
common curve as was already seen in numerical simulation of the random-Landau matrix model. 21 ) 
Furthermore, as a comparison, the behavior of the S(Q) dependence is also given for three families 
of P(s) functions based on Eq. (3): a) fixed (3 = 2 and variable a in the range [1, 2]; b) fixed a = 1 
and variable (3 in the range [0, 2]; and c) using (3 = 2^ and a = 7 + 1 with variable 7 in the range 
[0,1]. The first and the second families meat each other at the "semi-Poissonian" of the form 
P(s) = 2 ys 2 exp(— 3s) similar to that proposed in Ref 9. The family c is a Brody-like form that 
has already been seen to fit the evolution of the system in the S vs Q parameter space in Ref 21. 
As we can see our data indicate a similar behavior as the random-Landau matrix model. This is 
another nice example of the presence of one-parameter scaling since all data can be described by a 
one-parameter family of P(s) functions. This analysis suggests 7 0.54 around the critical point 
(E = 0). 

Right at the critical energy, however, the two parameter fit of P(s) to Eq. (3) is shown in Fig. 4. 
The fit was done for all available system sizes and resulted in the averaged values a = 1.56 ± 0.04 
and (3 = 1.99 ± 0.03. The errors were determined by the fitting and averaging procedures. Within 
those errors the result is almost the same as that obtained by Ono et al. for the Landau model in the 
lowest Landau band 14 ) (a = 1.53, j3 = 2.19). This contradicts the prediction that a = 1.21. 10 ' 11 ) 

In Ref 18 it was argued that the tail should show an exponetial decay. The results of a fit to a 
function exp(— ns) are shown in the inset of Fig. 4. We get k = 3.95 ± 0.6, which is comparable 
to the result by Batsch and Schweitzer 28 ) (n ss 4.1) and fits the prediction k = 1/(2%), with 
X = 0.124 ±0.06. 3 ) 

§4. Conclusion 

In conclusion we have presented the statistical analysis of spectra of pseudo-energies taken at 
different electron energies obtained in a network-model of the quantum-Hall effect. We could 
confirm that those pseudo-energy spectra show the same statistics as those of the energy eigenvalues 
of other models of the quantum Hall systems. This was done by investigating the shape of the level 
spacing distribution function P(s) and the scaling of its moments Jq = (s 2 ) /2, q = (s 2 )" 1 and 
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— (sins). Furthermore, although the overall shape of P(s) at the critical point has yet to be 
explained, we observe that level repulsion leads to the same quadratic behavior for small s as in 
the metallic regime, whereas the tail of the P(s) distribution falls off more slowly than in the GUE 
case. This is due to the multifractality of the critical wave functions and is also reflected in the 
finite compressibility (0 < x < 1). 
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Fig. 1. The P(s) spectrum for energies ranging from E — to E — 3. The smooth curves show the behavior for 
Poisson and GUE. The inset shows the same in double logarithmic plot. 



Fig. 2. The one-parameter dependence of the scaling variable Jo on L/£(E). The inset shows Jo as a function of 
energy E for different system sizes. At the critical point E c = all the curves meet at the same point. 



Fig. 3. The parametric dependence of S vs Q for our numerical results compared to three families of P(s) functions 
as described in the text. 



Fig. 4. The P(s) distribution for E — E c — 0. We can see the data for L = 20, 24, 30, 40, 50. All data points are 
positioned on the smooth curve that is a fit to Eq. (3), where a = 1.56 ±0.04 and /3 = 1.99 ±0.03. The inset shows 
a fit of exp(— ks) to the data for L — 40, k = 3.95 ± 0.06. 
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